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Abstract 

We study a quantum particle in a tilted honeycomb lattice in the tight-binding approximation. 
First we discuss the particle eigenstates, i.e., the stationary Wannier-Stark states. These states are 
proved to be extended states for the rational directions of the static field and localized states for 
the irrational directions. We find energy bands of the extended states and analyze the localized 
states. It is shown, in particular, that the localized 'honeycomb' Wannier-Stark states are chaotic 
states with irregular dependence of the localization length on the static field magnitude. Second 
we discuss Bloch oscillations of the quantum particle. Irregular Bloch oscillations for irrational 
directions are observed. 



I. INTRODUCTION 



Wannier-Stark states (WS-states) are eigenstates of quantum particle in a tilted lattice. 
Strictly speaking WS-states are resonances and have complex energies. However, for weak 
static fields they can be approximated by the stationary states with real energies (the single- 
band or tight-binding approximations). In the past two decades WS-states and related 
problems of Bloch oscillations (BO) and interband Landau-Zenner tunneling (LZ-tunneling) 
were readdressed in a number of fascinating laboratory experiments with cold atoms in 
(quasi) ID optical lattices and the light in ID arrays of optical waveguides, see Refs. [l-8| 
to cite a few of dozens relevant papers. These experiments stimulated theoretical studies 
that resuj^ed in essential progress in the theory of WS-states in one-dimentional systems, 
see Ref . 



ill for a review. 



An interesting extension of the theory of ID WS-states refers to 2D tilted lattices 10-121. 
It was argued in Ref. [10[ and confirmed later on in the experiment [7]] that WS-states in 
a 2D lattice are sensitive to the direction of the static field relative to primary axes of the 
lattice. Unfortunately, for square lattices considered in the above cited papers this effect is 
well pronounced only in the strong-field regime, where metastable nature of WS-states plays 
major role. In the present work we analyze WS-states and BO for a quantum particle in a 
honeycomb lattice. We will show that for this lattice the non-analitic angular dependence 
of WS-states is seen already in the weak-field regime, where the metastable WS-states can 
be approximated by the stationary states. This feature of the stationary 'honeycomb' WS- 
states has direct consequences for Bloch dynamics of the system that becomes qualitatively 
different for rational and irrational field directions defined later on in the text. It should 
be mentioned that BO in the honeycomb lattice were addressed earlier with respect to 
conductivity of graphene nano-ribbons [l^ [l^. However, in these works the electric field 
was aligned with the ribbon axis and, thus, the alignment effects were not discussed. 

This paper consists of two parts devoted to analytical and numerical analysis of WS- 
states, Sec. ini and BO of a quantum particle in the honeycomb lattices. Sec. IIIII The main 
results of the work are summarized in Sec. IIVI 
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II. HONEYCOMB WANNIER-STARK STATE 



In the standard presentation with two sublattices the tight-binding Hamiltonian of a 
quantum particle in the honeycomb lattice reads 

where J is the hopping matrix element, R denote coordinates of A-sites, and are three 
vectors that point from a A-site to the nearest B-sites. If a static field is present, this 
Hamiltonian should be complimented with the Stark term, 

H = Ho + E(F, R)aUR + E(F, RO^'^'R' (2) 

R R' 

where F is the field vector and R and R' are coordinates of the A- and B-sites, respectively. 
We are interested in the eigenstates of the Hamiltonian ([2]), i.e., in the stationary honeycomb 
WS-states. 

We begin with recalling general result concerning the structure of WS-states in a two- 
dimensional lattice of arbitrary geometry: For the field F parallel to a vector pointing from 
one lattice site to any other site WS-states are extended states in the direction orthogonal to 
F. These field directions can be labeled by two co-prime numbers r and q and for this reason 
are turmed rational directions. For example, for the square lattice the rational directions 
are given by Fy/F^ = r /q or 

tan6' = r/q , (3) 
while for the honeycomb lattice these are 

tan^ = y3^^. (4) 
q + r ^ ' 

The spectrum of WS-states for a rational direction (r, q) consists of infinite number of equally 

spaced energy bands. A particular feature of the square lattice and other simple lattices 

(like, for example, a triangular lattice with the hexagonal symmetry [15|) is that these 

energy bands are fiat for almost all rational directions. Because of this feature there is no 

qualitative difference between rational and irrational directions for the square lattice in the 

tight-binding approximation. In what follows we show that the Wannier-Stark energy bands 

of the honeycomb lattice have finite widths already in the tight-binding approximation. 

This makes a crucial difference between the honeycomb lattice and the square or triangular 

lattices, as well as between rational and irrational directions for the honeycomb lattice. 



A. Rational field directions 



For the rational directions WS-states in the honeycomb lattice are labeled by the ladder 
number n, the transverse quasimomentum k, and the sublattice index i. Correspondently, 
the energy spectrum is given by 

(k) = + dFn + e^^ (k) , (5) 

where 

1 ~ 3F 1 

Vr^ _l_ 2a yr^ — rq + q"^ 

and e*^*^(/t) is a periodic function of k. We calculated the spectrum ([5]) by adopting the 

method of Ref. 17 -l^. In brief, we map the honeycomb lattice into a square lattice with 

two sublattices. For this square lattice the static field is characterized by the vector F and 

the rational directions are given by Fy/F^ = r/q. Next we introduce another square lattice 

of the period d, which includes the previous lattice as a sublattice, and rotate it to align its 

X axis with the vector F. Finally we use substitution where the wave function is a plane 

wave along the y axis. After this sequence of transformations we end up with the system of 

two coupled ID equations, 

- J(e-*''"''^f_, + e'^'-^^f.^ + e'("-'?)"'^V'f-._,) + (dFj + E^)^ljf = Eijjf , (7) 

where E^ = (2/3) (r + q)dF if one sets E^ = 0. 

We solved Eq. ((Tj) numerically for different F and (r, g). Without any lost of generality 
one can restrict 6 to the interval < < 7r/3 which means g > r > 0. Examples of the 
spectrum (|5]) are given in Fig. [1] where the chosen energy intervals include three band with 
i = 1 and three bands with i = 2. In addition to Fig. [H Fig. [2] shows the width of the energy 
bands, 

A = maxe(/«) — minef/s;) , (8) 

as the function of F for the field direction 6 = n/G (asterisk) and some other directions. 
[Since A is independent of i, we drop the sublattice index in Eq. ([H]).] From the depicted 
numerical data one draws the following conclusions: The bands are well separated only in 
the limit of large F; With decrease of F the band width A monotonically grows while the 
distance dF between bands monotonically decreases and, for some Fcr ~ J, the bands almost 




FIG. 1: Fragments of the spectrum ([5]) for (r, 
field magnitude F. 
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FIG. 3: Widths of the energy bands for different rational directions for F = 2. 

touch each other. At this critical field magnitude A takes its maximal value; After reaching 
the maximum A shows monotonic decrease, where bands become flatter; This decrease is 
followed by erratic oscillations of the band width for small F. We note that in this region 
of small F the quantity ([8]) is not sufficient to characterize the spectrum because of rather 
complicated band pattern with many avoided crossings. 

We also studied the asymptotic behavior of A for F — )■ oo. Our numerical analysis reveals 
the dependence 

^-^.^ (9) 
where v increases with the increase of (r, q). For example, v = 1 for (r, g) = (1, 1), z/ = 2 for 
(r, g) = (1,2), etc. This result is similar to that for the band widths of the Landau-Stark 
states (eigenstates of a quantum particle in the Hall configuration) in the square lattice 
{igl - We believe that the power v = y{r,q) can be calculated analj^ically by adopting the 



perturbative approach of Ref. |19 |. 



Dependence of the band width A on the field direction 9 is depicted in Fig. [3l which 
presents the energy spectrum of the system (jS]) in the form of Hofstadter's butterfly. Namely, 
the figure shows energy bands for angles (jl]) with l<r<g<21, where the blue (dark 
grey) and green (light grey) colors correspond to i = 1 and « = 2, respectively. It is seen in 



Fig. |3] that the band widths progressively decrease with increase of (r, q) . Together with the 
estimate iQ this means that for irrational directions the spectrum is pure point and, hence, 
WS-states are localized states. 

B. Irrational field directions 

First we check that WS-states for irrational directions are localized states. In Fig. H] we 
compare two eigenstates of the Hamiltonian ([2]) with nearly the same energy for a 'rational' 
6 = 71 / 2 ^ 1.57 (equivalent to 9 = 7r/6) and 'irrational' 6 = 3 — 7t/2 ^ 1.43. The figure 
shows the integrated probabilities Py = /|\I'(R)pdx, dashed line, and = / |\I'(R)pd|/, 
solid line. (From now on we do not distinguish A- and B-sites.) For rational 6 the WS-state 
is seen to be an extended state in the direction orthogonal to F, while for irrational 6 it is 
localized in both directions. 

An important characteristic of the localized WS-states is their participation ratio, 

(Ei^wr) ' ^ (10) 

which indicates how many sites are occupied by a given state. We note that the participation 
ratio (|TOl) is the same for any WS-state because different WS-states are related to each other 
by translations. In the other words, one can obtain the complete set of WS-states from one 
state (or few states, if the lattice consists of sublattices) by translating it across the lattice. 




60 120 60 120 
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FIG. 4: Examples of extended (left panel, 9 = vr/2) and localized (right panel, 9 = 3 — 7r/2) 
WS-states with nearly the same energies. The dashed and solid lines show integrated probabilities 
along and across the field, respectively. The field magnitude F = 1. 
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FIG. 5: Localization length of WS-state for irrational direction 9 = 3 — tt/2. The thick and thin 
lines show the square root of the participation ratio (llOp of WS-states for the square and honeycomb 
lattices, respectively. The inset zooms into the region 1/F ~ 1. 

A remarkable feature of the honeycomb WS-states is that their participation ratio wildly 
oscillates if F is varied. The physics behind these oscillations is the following. Similar to the 
case considered in the previous subsection the (now discrete) spectrum of WS-state consists 
of two subsets that can be labeled by the subband index i or, what is the same, by the letters 
A and B. Correspondent ly, we have two families of WS-states 20|]. When F is varied the 



energy levels of A- and B-states non-monotonically move on the energy axis. If two levels 
of different symmetry come close to each other they develop an avoided crossing where the 
A- and B-states hybridize. As a consequence of the hybridization the function P = P{F) 
shows a local maximum, see inset in Fig. [5l It should be mentioned that to resolve all local 
peaks of P{F) (i.e, all avoided crossings in the spectrum) the step over F or, more precisely, 
over z = 1/F should be infinitesimally small. Figure [5] shows the function P = P{F) for 
a moderate step where only large peaks are resolved. Erratic oscillations with increasing 
density of peaks are clearly seen. This figure also reveals the expected average growth of 
the participation ratio when F is decreased. 

It is interesting to compare the localized honeycomb WS-states against the analytic results 
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for the simple square lattice, 



^(R) = Ji-n{J/F^)J^.u{J/Fy) , R = (/,m) 



(11) 



(here JJ^z) is the Bessel functions of the first kind). The square lattice has only one family of 
WS-states and, correspondently, the energy levels show no avoided crossings. Participation 
ration of the states ffTTj) is depicted in Fig. |5] by the thick line. Comparing two curves we 
conclude that WS-state in the honeycomb lattice are irregular or chaotic states, which are 
known to be sensitive to variations of the system parameters. The statistical analysis of 
these states in spirit of the Random Matrix Theory will be presented elsewhere. 

III. BLOCK DYNAMICS 

To study Bloch dynamics of the system ([2]) it is convenient to use the interaction repre- 
sentation with respect to the Stark term. This results in the time- dependent Hamiltonian 



where cuj = (F,rj) are the Bloch frequencies. Notice that the Hamiltonian f|T2l) commutes 
with the translation operator. Thus, when considering translation-invariant solutions of the 
Schrodinger equation, we can impose periodic boundary conditions. 

A. Delocalized initial state 

Consider initial state of the system \1/(R, t = 0) given by an eigenstate of the Hamiltonian 
([T]), i.e., by the Bloch wave \l/k(R) with the quasimomentum k. Then the time evolution of 
this state is naively expected to obey the equation \1/(R, t) ~ \l/k+Ft(R), that is known as 
the acceleration theorem. This simple Bloch dynamics, however, is complicated by the fact 
that the spectrum of consist of two subbands. 



with the energy gap vanishing at the Dirac points. Thus one generally can not avoid inter- 
band LZ-tunneling and the actual time evolution is given by the equation 

vl/(R,t)=f]Q(t)vl/«p,(R), (14) 




(12) 




(13) 
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FIG. 6: Populations of two Bloch subbands as the functions of time, upper panel, and the energies 
()13p for k' = k + Ft, lower panel. Parameters are 9 = tt — 3 and F = 1. 

where i is the Bloch subband index. The upper panel in Fig. [6] shows a typical dynamics of 
the coefficients Cj(t) and the lower panel depicts the energies f lT3|) at k' = k + Ft. It is seen 
that LZ-tunneling predominantly takes place when k' appears in the vicinity of the Dirac 
points. 

We would like to mention that the discussed Landau-Zener transitions may completely 
smear the periodic or quasiperiodic dynamics of quantum observables that is usually asso- 
ciated with BO. Figure [7] compares dynamics of p^, the projection of the mean momentum 
on the X axis, for F = 1 and = (upper panel) and 6 = n — 3 (lower panel). In the former 
case the particle trajectory in the (quasi)momentum space goes between Dirac points and 
LZ-tunneling can be neglected for chosen F. Thus the upper subband remains unpopulated 
and we observe nice periodic oscillations of px- In the latter case the upper subband gets 
populated independent of how small F is. As a consequence, px shows irregular oscillations 
where one hardly recognize the former periodic BO. 
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FIG. 7: Projection of the mean momentum on the x axis as the function of time for F = 1 and 
= 0, upper panel, and 6* = vr — 3, lower panel. 

B. Localized initial state 



Wave-packet dynamics in a 2D lattice with two Bloch subbands was considered earlier 



in Ref. 
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12| . The wave packet was found to have a tendency to spread in the direction 



orthogonal to F, while in the direction parallel to F it shows oscillatory dynamics. In this 
work we consider the limiting case of a localized initial packet where only one site is populated 
at t = 0. Also we will discuss dynamics in terms of WS-states instead of discussing it in 
terms of Bloch states as it was done in the above cited papers. To describe the wave-packet 
dynamics we introduce the time-dependent analogue of Eq. (|TOl) . 



(15) 



R 



According to the results of Sec. [TTlwe expect qualitatively different dynamics of the partici- 
pation ratio (fTSl) for rational and irrational field directions. 

The dashed line in Fig. |H] shows P(t) for F = 1 and rational direction ^ = 0. The 
participation ratio exhibits oscillatory dynamics superimposed with a linear increase in the 
mean value. This linear increase is due to ballistic spreading of the packet in the direction 
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FIG. 8: Participation ratio (jlSp for 9 = (dashed line) and 9 = tt — 3 (solid line). The initial 
condition corresponds to population of a single site in the center of the lattice. The value of the 
static force F = 1. 

orthogonal to F. The rate of baUistic spreading is obviously defined by the width A of 
the Wannier-Stark bands, while the characteristic frequency of oscillations is given by the 
distance between neighboring Wannier-Stark bands. 

The case of irrational direction = tt — 3 is depicted by the solid line in Fig. |8l Now P{t) 
saturates at some level defined by the characteristic localization length of the honeycomb 
WS-states. The wave-packet simulations also confirm l/F"^ scaling law for saturation level 
of P{t) that follows from the scaling law for the participation ratio of the localized 

honeycomb WS-states. 

We would like to mention that for some choices of the system parameters we observed 
rather exotic wave-packet dynamics. One example is 6 = 7i/6 and F = 2.3. As it is 
seen from Fig. [H^b), for these parameters the energy bands are almost straight lines, that 
implies nondispersive dynamics. Indeed, in our simulations we observed the formation of 
two soliton-like packets that propagate with constant velocities in opposite directions. 
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IV. CONCLUSIONS 



We found WS-states of a quantum particle in a tilted honeycomb lattice and compared 
them with WS-states in the tilted square lattice that are known analytically. The comparison 
is done for both rational and irrational directions of the field vector F. 

For rational directions of the field defined in Eqs. ([3lll]) the energies of WS-states form 
energy bands. For the square lattice these bands have zero width, excluding the case where 
F is aligned with one of two primary axes. This prohibits any transport in the system if the 
vector F is misaligned with a primary axis. Unlike this situation, for the honeycomb lattice 
the bands have finite width for any rational direction. If the condition (jl]) is fulfilled, an 
initially localized wave packet spreads in the direction orthogonal to F. It is interesting to 
notice that this type of wave-packet dynamics exist in the square lattice only in the presence 
of a gauge field normal to the lattice plane [19 1. 

For irrational directions of the field the energy spectrum of WS-states is discrete and, 
hence, they are localized states. We found the localization length of the honeycomb WS- 
states, which we define as the square root of the participation ratio ([TU]) . to grow in average 
as 1/-F, i.e., in the same way as for the square lattice. However, in a smaller scale the 
localization length shows large fluctuations that do not present in the case of square lattice. 
This observation motivates us to put forward a hypothesis about irregular (chaotic) nature 
of the honeycomb WS-states. This hypothesis is further supported by the irregular character 
of BO for irrational directions of the static field. 

The results of this work can be verified in laboratory experiments with cold atoms in 



honeycomb optical lattices 
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22|, honeyconib photonic crystals 23|, |2J], and microwave 



The common feature of these systems 



billiards with honeycomb array of scatterers 
is that they offer direct visualization of the 'wave function'. Another direction is ballistic 
conductivity of the graphene sheets. It is expected that the reported non-analytic angu- 
lar dependence of WS-states could strongly affect the conductivity. We reserve the latter 
problem for future studies. 
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